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Abstract. In this note we introduce three problems related to the topic of 
finite Hausdorff moments. Generally speaking, given the first n+1 (n £ NU{0}) 
moments, »o, ot\,..., ct n , of a real-valued continuously differentiablc function / 
defined on [0, 1], what can be said about the size of the image of 4^? We make 
the questions more precise and we give answers in the cases of three or fewer 
moments and in some cases for four moments. In the general situation of n + 1 
moments, we show that the range of the derivative should contain the convex 
hull of a set of n numbers calculated in terms of the Bernstein polynomials, 
x (1 — x) n+1 ~ k , k = 1,2, ...,n, which turn out to involve expressions just in 
terms of the given moments cti, i = 0,1,2, ...n. In the end we make some 
conjectures about what may be true in terms of the sharpness of the interval 
range mentioned before. 



1. Introduction 

We are studying here a problem from real analysis which can be roughly stated 
in the following way: 

given a continuously differentiable function whose first n moments 
are prescribed, what can be said about the image of the derivative of 
this function ? 

One of the tools that we will use is the following classical so called first mean 
value theorem for integrals (see Section 30.9 in p]). 

Theorem 1. Let h be a continuous function on [a, b] and g a non-negative Riemann 
integrable function. Then there exists a value c £ (a, b) such that 

rb rb 

h(x)g(x)dx = h(c) / g(x)dx. 



Moreover, if h{x) > h{c) (or h{x) < h{c)) for all x £ [a,b], then h(x) = h(c) for 
every x point of continuity of g and g(x) > 0. 

To introduce our hypothesis we let n £ NU{0} and let / be a continuously differ- 
entiablc function which satisfies the following Hausdorff moment type interpolation 
conditions: 



/ x k f(x)dx = a k , k = 0,1,2, ...,n, 
Jo 



a k e M. (1) 



Let us observe that given arbitrary moments a k the system ([lj leads to a linear 
one if / is a polynomial function. The main matrix of the resulting system is a 
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Hilbert matrix. This type of matrix is well know (see [2] , for instance) and has a 
non-zero determinant. 

Our investigation was motivated by a proposed problem in the College Mathe- 
matics Journal (|8 ) which requires one to show that if n = 2 and otk = k + 1, there 
exist ci,C2 G [0, 1] such that f'(c\) — —24 and /'(C2) = 60. It turns out that this 
problem was inspired by a problem of C. Lupu (see [6]) which referred to only two 
moments, ao = u\ = 1, and asked for a point c where f'{c) = 6. We wondered 
if these numbers were, in a certain sense which will be defined next, sharp. We 
will show that this is indeed the case in the next section (Theorem [5]) . Similar 
optimization questions, given the first n Hausdorff moments on [0, 1] or [—1,1], are 
customary subjects in the literature (see [5], [7]) We are going to formulate the 
following very general questions that are our main interest in this paper. 

Problem 1. For a fixed n and ctk as before, what is the largest range [A, B] such 
that [A, B] C Range(f') for every f a continuously differ entiable function on [0, 1] 
satisfying fifty? 

Problem 2. For a fixed n and as before, what is the biggest number L such 
that for every f a continuously differentiable function on [0, 1] satisfying fifty there 
exists some interval [a, b] with b — a = L that satisfies [a, b] C Range(f') ? 

We observe that in order to prove that [.A, 2?] 1S the answer for Problem [TJ it is 
necessary to show that [A, B] C Range(f') for every / a continuously differentiable 
function on [0, 1] satisfying (JT]) and that for every e > there exists fi and f r 
continuously differentiable functions on [0, 1] satisfying (JlJ and 

Range(fl) C (A — e, 00) and Range{f' r ) C (—00, B + e). (2) 

It is clear that if A and B give the answer in Problem [TJ then in trying to answer 
Problem [5] we must have L > B — A. If for every e > 0, one can find a function 
(/; = fr) that satisfies both conditions in ([2]), then the answer to Problem [2] is 
simply L = B — A. 

Another related problem here is to characterize the case L > B — A and calculate 
L in this case in terms of the afc's. Perhaps Problem 2 may be easier if one 
restricts the class of functions in consideration to something more manageable like 
polynomials of a certain degree. 

If we want to make the range of the derivative as small as possible, we just have 
to take moments that satisfy the necessary and sufficient condition for having a 
solution to the system that results from having a linear function, say f(x) = u + vx, 
x e [0, 1], satisfying (TTJ): 

= a k , k = 0,1,2, ....n. 



k+1 k+2 
This is equivalent to 



rank 



1 h "o 

1 f 

? I 

3 4 



0:1 
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On the other hand, if we want to make the range of /' as big as possible, it 
makes sense to restrict our moments to a finite range, say [—1, 1]. We observe that 
the problem is homogeneous under dilations, so let us formulate a third problem. 

Problem 3. For a fixed n, what is the maximum of B — A such that [A, B] C 
Range(f') for every f a continuously differ entiable function on [0, 1] satisfying (QJ), 
the maximum being taken over all possible moments a& G [—1,1]? 

We will show in Section 2 that the answer to Problem [3] is 156 if n = 2, for 
the moments ao — 1, a\ = — 1 and a 2 = 1- We observe that if the answer to 
Problem [2] is zero, then the answer to Problem [3] is also zero. As suggested by one 
of the referees of our paper, one can ask similar questions about the range of /" or 
higher derivatives, assuming these exist. We will make some remarks about these 
questions and see how the results for the first derivative could be applied for higher 
derivatives. 



2. Small values of n 

We have a few complete answers to Problem Q] for small values of n (n < 3). 
First, let us study what happens with n = 0. If we take g{x) = 1 — x and h = f 
in Theorem [TJ using integration by parts, we get 

f(c,)l = / f'(x)(l - x)dx = -/(0) - f 1 f(x)(-l)dx = a Q - /(0), Cl G (0, 1), 
1 Jo Jo 

or 

/'(ci)= 2(ao-/(0)), Cl G(0,l). 
If /(0) = a, then we can take f{x) = a+ (2ao — 2a)x and observe that in case n = 0, 
there exists a function such that J f(x)dx = ao and Range(f') = {2ao — 2/(0)}. 
This gives us the following simple answers to Problem [T] and Problem [2j 

Proposition 1. For n = 0, there is no A and B that satisfy the requirements of 
Problem]]] The answer for Problem^ (n = 0) is L = 0. 

Let us continue the analysis in the case n — 1. We can apply Theorem [T] to 
g(x) = x(l — x) and h = f , x G [0, 1]. Then, a similar calculation gives that for 
some C2 G (0, 1), 

1 f 1 

f(c 2 )- = - / f(x)(l - 2x)dx = 2a! - a , => /'(<*) = 6(2 ai - a ). 
If we apply Theorem [T] to h = f and g(x) = (1 — x) 2 instead, 

f'(c 3 )\ = f f'(x)(l-x) 2 dx = -/(0)- / f(x){2x-2)dx = 2(o -ai)-/(0), c 3 G (0, 1), 
Jo Jo 

or 

/'(c 3 ) = 6(a - ai) - 3/(0), for some c 3 G (0, 1). 

So, if we take a — 2(2ao — 3ai) and f{x) = a + mx where m = 6(2ai — ao) = 
6(ao — ai) — 3a, we get a function which will give us what we need in this case, 
and therefore provide a similar answers to our problems. 
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Proposition 2. For n = 1, we can take A = B = 12ai — 6ao to satisfy the 
requirements of Problem^ The answer for Problem^ (n = 1 ) is L = 0. 

The case n = 2 is getting a little more interesting; it is essentially non-trivial 
and at the same time pretty surprising. We have a definite answer to Problem [T] 
and Problem [3] and we show some inequality for L in Problem [5] 

Theorem 2. For n = 2, if Aq := 602 — 6«i + ao > 0, the values 

A := 12(4ai - a - 3a 2 ) and B := 12(3a 2 - 2a x ) 
satisfy the requirements of Problem [7] and if A < then one needs to switch 
the values of A and B above in order to solve Problem [7J // Ao = 0, the values 
A = B = 12(3a2 — 2ai) answer Problem^]] and L = answers Problem^ 

Proof. First, let us show that A and B are always in the range of the derivative. 
This is done as we have seen before by setting in the Theorem [TJ h = /', and 
g{x) = x(l — x) 2 > (x G [0, 1]). Indeed, we have J* g(x)dx = -^ and 

1 ,1 
/' \x)g(x)dx — f(x)g(x)\l — / f(x)(l — 4x + 3x 2 )dx = 4ai — ao — 3a2. 



Jo 
Hence, for some C4 we must have /'(C4) f Q g(x)dx = 4ai — ao — 3a 2 which in turn 
gives /'(c 4 ) = 12(4ai — «o ~ 3^2) = A. Similarly, for g(x) = x 2 (l — x), (x 6 [0, 1]), 
one finds that f Q g(x)dx = still holds true and 

/ f'{x)g{x)dx= I f{x){2x — 3x 2 )dx = 3a2 — 2a±. 
Jo Jo 

This insures that B = 12(3a2 — 2ai) is also in the range of /'. Because /' is 

assumed to be continuous we get that the whole interval [A, B] or [B, A] is contained 

in the range of /'. 

From here on, we are going to work under the first assumption (Ao > 0) which 
is equivalent to A < B (B — A = 12Ao). To show that A and B are sharp bounds 
we begin with B by constructing a spline function St for t € (0, 1), defined by 



a + bx + cx for x € [0, t] 



st(x) = 



m + nx if x G [t, 1], 

where a, 6, c, m and n are determined by the conditions L St{x)dx — ao, Jq xst{x)dx 

ai, J Q x 2 St (x) dx = a2 and the restrictions necessary to insure that st is continu- 
ously differentiable at x = t. 

In order to add an intuition element we included here the graphs of Si/10 and 
its derivative for ao = 1, a.x = 1 and a2 = 2. 

It is easy to see that s t is continuously differentiable at x = t if and only if 
m = a — ct 2 and n = b + 2ct. This gives the new expression of St just in terms of 
a, b and c: 



a + bx + cx 2 for x G [0, t] 



st(x) = { (3) 
(a - ct 2 ) + {b + 2ct)x if x e [t, 1]. 
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(a) s_ 



04 06 



(b) s\ 
TO 



Figure 1 . Case n = 2 and oik = k + 1 



One can check that the moment restrictions reduce to the following 3x3 relatively 
simple linear system of equations in a, b and c: 



a b 
<2 + 3 



12 



t + 1 I c = ao 



b 
4 



30 



3 ~ 2 ' ' 



(4) 



a 2 . 



2 3 

It is clear that we have a unique solution for this system at least for infinitely many 
values of t since the main determinant of the system is a polynomial in t of degree 
at most 5. Let us observe that 



'b + 2cx for x e [0,t] 



ax 



(5) 



b + 2ct if x G [t, 1]. 

We observe that if c > 0, the maximum of this function is b + 2d. With a little 
work one solves the system (Q| and finds that 



30A , „ 12(6t 2 ai - 3t 2 a - 20«i + 30a 2 + 5to - 15ta 2 ) 

c = -=-; — ^ r : o + 2ct = . 

t 3 {6t 2 - 15* + 10) 6t 2 - 15i + 10 

It is clear from these expressions that under our hypothesis c > for every t > 
and that B — 12(3a 2 — 2ai) = lim(6 + 2ct) which proves that B is sharp. In a 

similar way one can prove that A is sharp by taking a spline s t which is first a linear 
piece on [0,*] and a quadratic piece on [t, 1]. It turns out that the calculations are 
very much similar to the ones above with the only difference that this time we let 
t approach 1. 
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However, we are going to show that the lower bound A is sharp by using an 
invariance principle here by doing a "change of variable" so to speak and considering 
how the Problem [T] changes from / to g where g(x) = f(l — x), x G [0, 1]. Let us 
denote by ctk(f) the k th moment for the function /. One can see that 

a (g) = a (f), ai(g) = a (f) - ax(f) and a 2 (g) = a (f) - 2ai(/) + o 2 (/), 

and of course, the relations are symmetric with respect to interchanging / and g , 
i.e. 



«o(/) = "0(5), anif) = a (g) - 0*1(3) an d a 2 (f) = a (g) - 2a 1 (g) + a 2 (g). 

Let us observe that the hypothesis that Ao(/) > is in fact invariant under this 
change: 

6o 2 (/) - 6«i(/) + a (f) = 602(5) - 601(5) + a (g) > 0. 
By the first part of our proof, we see that 
B{g) = 12(302(5) - 201(5)) = 12[3(o (/) - 2oi(/) + o 2 (/)) - 2(o (/) - oi(/)), 
or 

B(5) = 12(3o 2 (/)-4oi(/)+a (/), 

is a sharp bound for the range of 5'. Since g'{x) = — f'(l — x) we see that the 
range of 5 is just the range of / reflected into the origin and vice versa. Hence, 
A(f) = -B{g) = 12[4oi(/) - 3o 2 (/) - o (/)] is a sharp lower bound for /. The 
rest of the statements of the theorem follow from what we have shown so far. ■ 



Corollary 1. In the case n = 2, in respect to Problem^ we have 

12|A | < L < 32|A |. 
The maximum required in Problem\3\is 156. 

Proof. The first part is a simple consequence of the fact B — A = 12Ao and the last 
part follows from the fact that |A | = 6o 2 — 6oi+o | < 13 if o, G [—1, 1], i — 0, 1,2. 
To show the inequality L < 32|Ao| we employ the same idea by constructing spline 
which is symmetric around 1/2: 



'a + bx + c{x - 2tx - t 2 + t - 1/4) for x S [0, 5 - 1] 
st(x) = < a + bx + ex 2 if x 6 [5 — t, | + 1] (6) 

a + bx + c(x + 2tx - t 2 - t - 1/4). 

This spline is continuously diffcrcntiable on [0, 1] and depends on three parameters 
which if determined from the constraints given by the moments we get 



120A n 



i(15-40i + 48i 2 ) 



> 0,for all t G [0,11, 
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and which shows that the minimum and the maximum of the derivative of §t is 
attained on the linear pieces. One can see that the difference between these two 
values is actually Act and so letting t — >• we get that 



L < lim Act 



32|Ao| 



One can use the same techniques to show that for three moments, assuming 
the second derivative exists, the range of the second derivative should contain 3(Wo 
and this is sharp because a polynomial of degree two exists solving the moments 
problem. 

The case n = 3 is even more interesting and a lot more complicated. First of all 
we have at least three new possible values that we need to add to the range of /': 



C = 20(4a 3 -3a 2 ), D = 60(3a 2 -2a 3 -ai), and E = 20(4a 3 -9a 2 + 6ai -a ) (7) 

obtained from Bernstein polynomials, gi(x) — x 3 (l — x), 32(2;) = x 2 (l — x) 2 and 
gs(x) = x(l — x) 3 respectively. 

Proposition 3. Given A and B as defined in Theorem^ we have the inclusion 

[min(A,B),max(A,B)] C [min(D,C,E),max(D,C,E)]. 

Proof. Let us observe that g\{x) + g%{x) — x 2 (l — x) and g%{x) + 33 (x) = x(l — x) 2 . 
Differentiating and integrating against f(x) we get the relations j^D + = j^A 
or A — l-D+l-E 1 and similarly B = |C+|D. These two convex linear combinations 
are enough to conclude the desired statement. I 
Of course, this proposition can be generalized to an arbitrary n. So we expect 
that the interval that answers Problem [T] contains the convex hull of the numbers 
constructed as usual, i.e. 



D, 



tif{x)±{x k {l-x)^- k )dx 



k = 1,2, 



(8) 



/„ x k (l - x) n + 1 - k dx 

given by the highest degree Bernstein basis polynomials possible. 
We observe that if we define Ai := 10a 3 — 12a2 + 3ai, then 

C = E + 20A and D = C - 20Ai. 

Hence, we observe that if we have Ao > and Ai < for instance, then D > C > E. 
Therefore, in light of Proposition [3l the candidates for the two values needed to 
answer Problem Q] are A = E and B = D under the given assumption. In fact, for 
various other situations we believe that the values A and B that answer Problem [1] 
are given for each case in the following table 



No 


Hypothesis 


A 


B 


(i) 


A > 0, Ai < 


20(4a 3 — 9ai2 + 6ai — ao) 


60(3a2 — 2a 3 — ai) 


(ii) 


< A < Ai 


60(3a2 — 2a 3 — a±) 


20(4a 3 - 3a 2 ) 


(iii) 


< Ai < A 


20(4a 3 — 9a 2 + 6ai — ao) 


20(4a 3 - 3a 2 ) 


(iv) 


A < 0, Ai > 


60(3a2 — 2a 3 — ai) 


20(4a 3 — 9a2 + 6ai — ao) 


(v) 


Ai < A < 


20(4a 3 - 3a 2 ) 


60(3a2 — 2a 3 — ai) 


(vi) 


A < Ai < 


20(4a 3 - 3a 2 ) 


20(4a 3 — 9a2 + 6ai — ao) 



8 



EUGEN J. IONASCU AND RICHARD STEPHENS 



where A = 72 



Oil 
«2 



and Ai = 720 



Q 2 



We have the following partial result along these lines. 

Theorem 3. For n = 3, the upper bound of (i) and the lower bound of (ii), in the 
table above, are correct. If Ao = Ai = then A = B = 20(4a3 — 3a?2) and L = 
solves Problem^ and Problem^ 

Proof. First of all let us observe that the cases (iv) , (v) and (vi) follow from (i) , (ii) 
and (iii) respectively by simply changing / into — /. This simple transformation 
changes basically the order of A and B. It is easy to see that Ao = Ai = implies 
the existence of a linear map that has the given moments and so A — B and L = 0. 
Hence in what follows we will assume that Ao ^ or Ai / 0. 

Based on the invariance principle that we used in the proof of Theorem[2]we need 
to show the sharpness of only the upper bound in (i). Indeed we observe that if 
g(x) = f(l—x), x 6 [0, 1] then one can check that the hypothesis Ai(/) < changes 
into Aq(</) < Ai(<?). Also, the hypothesis < Ai < A , is actually invariant under 
this change. One also needs to take into account that the bound D is invariant 
under this transformation but C and E interchange: 

D(f) = D(g), C(f) = E(g), and E(f) = C(g). 

So, let us begin with case (i) and show that B = 60(3a2 — 2a 3 — at±) is sharp. 
For every t G (0, 1/2), consider a spline function si it which is quadratic on [0,t], 
linear on [t, 1 — t] and another quadratic on [1 — t, 1]. The constraints of having 
this spline a continuous and differentiable function give us a similar form for si.t 
to the one constructed in the proof of Theorem [2] in equality ([3]), in terms of four 
free parameters a, 6, c and d: 



bx 



cx' 



for x G [0, t] 



Si.tO) = { (« - ct2 ) + {b + 2ct)x if x G [t, 1 - t] 

d(l -t) 2 + a - ct 2 + [b + 2ct - 2d(l - t)]x + dx 2 if x G [1 - t, 1]. 

The four parameters are then determined by imposing the four linear constraints 
given by the moments. The resulting system is 



t 3 

+ ^ + C-7r-t 2 +t)c+—d- 



3 



a 



i 4 


2t 


i 2N 


12 


+ 3 




t 


t 2 


~) 


2 " 


-3 + 


30 J 


2t 


t 2 




~5 


4 " 


f 60, 



d = otx 



(9) 



+ — }d 

3 6 30, 



a 2 



*_ 
10 



60 



a 3 . 
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As we have observed before, the system has a unique solution for infinitely many 
values of t 6 (0, 1/2), since the main determinant of the system is a polynomial in 
t of degree at most 11. Because the derivative of s± t t is given by 

'b{t) + 2c(t)x for x e [0,i] 

s i,ti x ) = { b (t) + 2c (t)t if xe[t,l- 1] 

b(t)+2c(t)t-2d(t)(l-t-x) if xG[l-t, 1]. 
One can use a symbolic calculator and check that 

lim b(t) + 2c{t)t = 60(3a 2 - 2a 3 - at), 

which is one necessary fact to prove the sharpness of B. Also, we need to check 
that for most of the values of t, b(t) + 2c(t)t is a maximum of the derivative of Si^. 
For this end, it is enough to check that c(t) > and d(t) < for small values of 
t. Again, one can compute lim^o c(t)t 3 = Aq — A\ > under our assumption in 
case (i) (unless both numbers Aq, Ai are zero). Also, the limit of d(t)t 3 as t — > 
turns out to be equal to 3Ai < 0. If Ai = we know that Ao > 0. In this case we 
have lim t ^o d{t)t 2 = -|A < 0. 

Using the duality via g(x) = f(l — x), we see that A = 60(3a2 — 2a 3 — a%) is a 
sharp lower bound in the case (ii). I 

In the case n = 3, assuming the table before Theorem[3]is correct, with respect to 
ProblemEl we have either L > 20|A |, L > 20|Ai|, or L > 20|A - Ai|, depending 
upon the hypothesis in which the moments fall into as classified in Theorem[3] The 
maximum required in Problem [3] is 760 which is attained for ao = 1, ot\ — — 1, 
Q2 = 1, and 03 = — 1. We wonder if alternating the signs of the moments and 
setting them ctk = (— l) fc will always give the maximum in Problem[3l 

For higher derivatives we can show that U := 120(3ai — 12a 2 + lOa^) and 
V := 120(ao — 9ai + 18a2 — 10«3) are in the range of the second derivative. It does 
not seem to follow from our Theorem [2] applied to /' that these values are sharp, 
although the same idea of using a spline formed by a a cubic and a quadratic may 
work. 



3. Higher values of n 

We have noted the following statement after the proof of Proposition [3] 

Theorem 4. Given a continuously differentiable function satisfying the Hausdorff 
moments constraints {T]) (n > 2), then the range of the derivative contains the 
interval [A n ,B n ], where A n — minjDfclfc = 1,2, ...,n} and B n = max{£>fc|£; = 
1, 2, n}, with Dk given by (0). Moreover, [An, B n ] C [A n +i, -B„+i] for all n > 2. 

Proof. The first part follows with the same technique we have employed over and 
over here using Theorem [T] For the second part we are observing that the integrals 
which appear in the denominators of ([8]) , are actually the well known values of the 
beta function, i.e. B(k + l,n + 2 — k) = J x k (l — x) n+1 ~ k dx. Using the established 
formula for B( , ), we see that 
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B(k + l,n + 2-k) 



T(k + l)r(n + 2-k) _ k\(n + 1 - k)\ 



T(n + 3) 



(n + 2)! (n + 2)("+ 1 )- 



This gives us a new expression of the D k ^ n which is basically in terms of the genuine 
Bernstein basis polynomials, i.e. b un = (^\x"(l — x) n ~", v — 0, 1, ...,n: 



rl d 

Dk, n ■= -(n + 2) / f(x) — (b kin+1 )dx, k = 1,2, ...,n. 



n 



(10) 



It is easy to check that x fe (l - x) n+1 k + x k+1 (l - x) n k = x k (l - x) n ~ k which 
basically gives the convex combination formula 



n + 2-k k+1 
D kn = — D k<n+ i H —D k+1 . n+ i, k = 1,2, ...,n, n > 1. 



n + 3 



n + 3 



These expressions imply the second claim of the theorem. ■ 
Let us observe that ([TT]) implies the following form for D kiH 

D k>n :=(n + l)(n + 2) / f(x){b k<n - b k - 1>n )dx, k = l,2,...,n. (11) 
Jo 

which provides a simple way of computing _Dj.,s in terms of the moments ao, cx\,..., 
a n . In what follows we will describe yet another way of doing these computations, 
and for that purpose we generalize first the definitions of Ao and Ai in the following 

way 



A fc = 



(fc + l)(fc + 2) 2 (fc + 3) 2 (fc + 4) 



i l 

fc+1 fc+2 



fc^-2 a k+l 
fc+3 k+4, a k+2 



, k = 0,1, 2, 



or simply 



(fc + 3)(fc + 4) n . n . (fc + l)(fc + 2) 
A fe = i '-a k+2 - (k + 2){k + 3)a k+1 + ± ^ J -a k , k>0. 

There are some relations between the D' k s and A' k s in general which we will include 
in the next proposition. 

Proposition 4. For k > and n > 2, we have in general 

Afc = 2(fcT3) ' fc "°- (12) 

71-2 

£> n ,„ = 6(2a 1 -a ) + 2^(fc + 3)A fe . (13) 

fc=0 

Moreover, with the definitions of A n and B n from Theorem^ A n — B n if and only 
if Ai — for all i = 0,1, 2, n — 2, if and only if there exists a linear function 
with moments ct k . 
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Proof. Let us observe that we can simply write 

x k = x k+1 +x k (l-x) = x k+1 + ^jh,k+i 

fca fc _ 1 = (fc + l)a fc - (fc + 1) 1 (fc + 2) D fc , fc . 
Using this last formula we can calculate the expression of A&: 



cA fc = 



1 

fe+i 
1 

fc+2 
1 

fc+3 



1 

fc+2 
1 

fc+3 
1 

fc+4 



f±f«fc+l 



(fc+l)(fc+2)(fc+3) 

afc+i 



fc+2 „,. I -Pfc+2.fc+2 
fe+3 «fc+l T ( fc+ 3)2( fc+4) 



1 


1 


fc+1 


fc+2 


1 


1 


fc+2 


fc+3 


1 


1 


fc+3 


fc+4 



-Pfc + l,fc + l 
' (fc+1) (fc+2) (fc+3) 



-Pfc + 2 fc + 2 

(fc+3) a (fc+4) 



where c 



(fc+i)(fc+2) :i (fc+3)^(fc+4) • This last identity implies the formula ()12|) . For 
the second part of our statement we observe that A, = for all i = 0,1, 2, ...n — 2 
if and only if there exists a linear function / with moments at- In this case the 
range of the derivative of / consists of only one point and therefore by Theorem [4] 
we must have A n = B n . For the converse, again using Theorem |4] we obtain that 
all Dij, 1 < i < j, 1 < j <n, have identical values and so by (fl"2"T) we get A& = 
for all k = 0,1,..., n-2. 

Finally, let us observe that the equalities in (fT2")l provide a telescopic sum for 
Dk,k which allows one to arrive at formula (fT3|). ■ 
The convexity relations can be used to calculate all the D' k s from the Dfc,fc and 
so formulae (|13|) provide a way of computing all the D' k s in terms of determinants 
A,. 

For the case = k + 1, k = 0, l,2,...,n, we calculate A n and B n in a more 
precise way. This generalizes the problem in [?] . 

Corollary 2. Let n G N, n > 2, be fixed and f be a continuously differentiable 
satisfying {7Jj with = k + 1, k = 0, 1,2, ...,n. Then, the values of A n and B n as 
defined in Theorem are 

A n = -n(n + l)(n + 2), B n = (n + l)(n + 2)(2n + 1). 

Proof. Using the formula for A& we get 



(fc + 3) 2 (fc + 4) 
2 



(fc + 2) 2 (fc + 3)+ (fc + 1) ' (fc + 2) =3fc + 7, fc>0. 



Then using formula (|13[) we obtain 

n-2 



A»,n = 18 + 2^(/s + 3)(3fc + 7) = (n+ l)(n + 2)(2n+ 1). 



fc=0 



Now we can use the convexity relations and compute _D„_i.„: 
1 



= -n(n + l)(n + 2). 
Next, if one calculates -D n -2,n, some surprise appears: 



^n-l.n = ^ (( n + 2 )- D n-l,n-l - 



D n -2,n = g (fa + 2)D n _ 2 ,n-l - (" - l)-D„-l, n ) = 0. 
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Because of the convexity relation, it is easy to see that all the other Dk. n , k < n — 2, 
are equal to zero. Therefore A n — — n(n + l)(n + 2) and B n = (n+l)(n+2)(2n+l). 
■ 

Putting together what we did so far we now can say that for oik = k + 1, the 
bounds above are sharp if n = 2 and the lower bound is sharp if n = 3. 

Theorem 5. For n > 2 fixed, with the definition of A n and B n as in Theorem^ 
if An < B n it is not possible to have L = B n — A n in Problem^ 

Proof. By way of contradiction let us assume that L = B n — A n . Hence, we can 
find a sequence of functions f m , continuously diffcrentiable, such that Range(f m ) C 
[A n — i, B n + i] and satisfying (|TJ). Since f' n can be considered in L 2 ([0, 1]) we can 
find a subsequence of f' m , say /' , weakly convergent to a function / € L 2 ([0, 1]). 
This implies that for every non- negative function g £ L 2 ([0, 1]), 



(A n - < I f' m {x)g{x)dx < (B n + -)\\g\\ 

m 



o 



where \\h\\i — L h(x)dx, h G i 1 ([0, 1]). Passing to the limit as — > oo, we get 



^n||fl||i< / f(x)g(x)dx < BnWgWi, S6l 2 ([0,l]),j>0. 







This implies that A n < f{x) < B n for a.e. x E [0, 1], by a standard measure theory 
argument. Since B n = for some k = 1,2, ...,n, and f m satisfy ([lj we can say 
that 

||s*(l-a:) n + 1 -*||i ||x fc (l-a;)™+ 1 - fc ||i 
Letting — > oo we obtain 

_ $lf{x)x k {l-x) n + x - k dx 



||a; fe (l 



\n+l— k 



1 



Hence, re- writing this yields 

[B n - f(x)]x k (l-x) n+1 - k dx = 0, 

which in turn implies, by what we have shown before about /, that f(x) — B n for 
a.e. x G [0, 1]. Similarly, we arrive at the conclusion f(x) = A n for a.e. x <E [0, 1]. 
Since we assumed A n < B n we clearly get a contradiction. Therefore, it remains 
that L > B n — A n . ■ 
This last theorem says that Problcm[T]and Problem[2]are completely different in 
nature. We must admit that we do not have a definite answer to Problem [2l other 
than the trivial case L = 0, in any of the particular situations we have considered. 
We leave that to the interested reader. 
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